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20.  ABSTRACT  (continued) 

\ 

'■'example  that  quantitatively  simulates  the  characteristics  of  a typical  tow-tank 
test  vehicle.  - 

The  maintenance  of  steady  laminar  flow  on  a heated  submerged  vehicle  requires  that 
the  external  surface  temperature  at  every  location  along  the  vehicle  be  equal  to  or 
greater  than  that  required  to  delay  transition.  The  ability  of  the  vehicle  to  re- 
cover laminar  flow  is  sensitive  to  the  amount  by  which  the  wall  temperature  exceeds 
the  minimum  required  for  laminar  flow.  If  there  were  no  temperature  excess,  the 
recovery  time  would  be  infinite.  When  there  is  excess  temperature,  the  recovery 
time  depends  on  the  extent  to  which  the  turbulent  event  has  lowered  the  skin  temper- 
ature, and  is  characterized  by  a time  constant  that  is  directly  proportional  to  the 
heat  capacity  per  unit  are  of  the  skin  and  the  heating  fluid  together  and  inversely 
proportional  to  the  laminar  heat  loss  rate  to  the  external  flow. 

For  the  calculated  cases,  if  an  encountered  patch  of  turbulence  extends  the  order 
of  three  or  more  vehicle  lengths,  the  recovery  time  for  reasonable  temperature 
excesses  is  of  the  order  of  one  to  three  characteristic  times  (1/2  to  2 minutes  for 
the  cases  calculated).  For  momentary  losses  of  laminar  flow,  recovery  occurs  very 
soon  after  the  turbulence  is  swept  past  the  trailing  edge. 
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ABSTRACT 


In  order  to  study  the  process  by  which  laminar  flow  is  re-established 
after  a heated  submerged  vehicle  encounters  a patch  of  turbulent  flow 
an  analytical  model  has  been  formulated  which  simulates  the  heat  trans- 
fer characteristics  of  the  external  skin  and  internal  heating  passage 
of  such  a vehicle.  This  model  considers  the  energy  equations  for  the 
external  skin  and  for  the  internal  heating  passage  and  the  coupling 
between  them  due  to  heat  transfer  to  the  external  flow.  Numerical 
solutions  to  this  coupled  set  of  equations  for  a number  of  turbulent 
catch  scenarios  have  been  obtained  by  an  explicit  finite-difference 
procedure  for  a two-dimensional  flat-plate  example  that  quantitatively 
simulates  the  characteristics  of  a typical  tow-tank  test  vehicle. 

The  maintenance  of  steady  laminar  flow  on  a heated  submerged  vehicle 
requires  that  the  external  surface  temperature  at  every  location  along 
the  vehicle  be  equal  to  or  greater  than  that  required  to  delay  transi- 
tion. The  ability  of  the  vehicle  to  recover  laminar  flow  is  sensitive 
to  the  amount  by  which  the  wall  temperature  exceeds  the  minimum  required 
for  laminar  flow.  If  there  were  no  temperature  excess,  the  recovery 
time  would  be  infinite.  When  there  is  excess  temperature,  the  recovery 
time  depends  on  the  extent  to  which  the  turbulent  event  has  lowered  the 
skin  temperature,  and  is  characterized  by  a time  constant  that  is  di- 
rectly proportional  to  the  heat  capacity  per  unit  area  of  the  skin  and 
the  heating  fluid  together  and  inversely  proportional  to  the  laminar 
heat  loss  rate  to  the  external  flow. 

For  the  calculated  cases,  if  an  encountered  patch  of  turbulence  extends 
the  order  of  three  or  more  vehicle  lengths,  the  recovery  time  for  rea- 
sonable temperature  excesses  is  of  the  order  of  one  to  three  character- 
istic times  (1/2  to  2 minutes  for  the  cases  calculated).  For  momentary 
losses  of  laminar  flow,  recovery  occurs  very  soon  after  the  turbulence 
is  swept  past  the  trailing  edge. 
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1.  INTRODUCTION 

In  recent  years,  there  has  been  a resurgence  of  interest  in  attempting 
to  improve  a vehicle's  performance  by  increasing  the  length  of  laminar 
flow  in  its  external  boundary  layer.  By  delaying  the  onset  of  turbu- 
lence (known  as  Laminar  Flow  Control  or  LFC),  the  total  drag  of  the  ve- 
hicle may  be  reduced  which,  in  turn,  would  allow  either  the  use  of  a 
smaller  propulsion  system  for  a given  mission  or  for  an  extended  range 
or  top  speed  for  a given  propulsion  system.  A number  of  techniques 
for  increasing  the  extent  of  laminar  flow  have  been  identified  such  as 
vehicle  shaping  (wherein  a favorable  pressure  gradient  along  the  vehicl 
delays  the  tendency  for  the  laminar  boundary  layer  to  undergo  transi- 
tion to  turbulence)  and  suction  of  the  boundary  layer  fluid  through 
the  vehicle  skin.  For  submersible  vehicles  (as  opposed  to  airships) 
the  technique  of  heating  the  vehicle  surface  also  has  been  found  to 
delay  transition  because  of  the  temperature-viscosity  behavior  of  water 
This  latter  technique  has  the  added  attraction  that  waste  heac  from 
the  vehicle's  propulsion  system  might  be  made  available  to  heat  the 
vehicle  surface  rather  than  simply  being  exhausted  overboard;  thus  an 
increased  laminar  flow  length  might  be  obtainable  at  little  or  no  added 
expenditure  of  power  (Reference  1).  Experimental  research  is  in  pro- 
gress to  determine  how  much  of  the  theoretical  performance  increase  !s 
actually  achievable.  Systems  studies  are  also  being  conducted  to  deter 
mine  if  sufficient  waste  heat  is  available  and,  if  not,  how  it  might  be 
supplemented  without  seriously  affecting  vehicle  design  goals. 

One  aspect  of  the  Laminar  Flow  Control  area  that  has  not,  as  yet,  re- 
ceived much  attention  has  to  do  with  the  recovery  of  laminar  flow  once 
it  has  been  partially  or  totally  lost.  For  example,  any  submerged  ve- 
hicle will  occasionally  encounter  isolated  regions  of  turbulence  (the 
oceanographic  equivalent  of  clear  air  turbulence)  or  the  laminar 
boundary  layer  flow  may  be  disrupted  by  the  adherence  of  a particle 
to  the  vehicle  surface  - thus  acting  as  a temporary  boundary  layer 
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"trip".  There  is  even  the  possibility  that  such  a particle  may  simply 
disrupt  the  boundary  layer  flow  and  create  a turbulent  "spot  which  is 
swept  along  in  the  boundary  layer  across  the  vehicle's  surface.  The 
loss  of  laminar  flow  is  especially  crucial  for  a heated  vehicle  since 
turbulent  heat  transfer  rates  are  so  much  greater  than  laminar  rates 
i that  the  vehicle  can,  in  effect,  lose  its  overheat.  Depending  on  the 

vehicle's  heating  system  and  overall  heat  capacity,  a considerable 
amount  of  time  may  be  required  before  it  can  heat  up  again  and  re- 
establish laminar  flow.  The  modeling  and  calculation  of  such  an  occur- 

* rence  is  the  subject  of  the  present  report. 

Because  this  is  the  first  time  (to  the  authors'  knowledge)  that  the 
laminar  flow  recovery  problem  has  been  analytically  addressed  (see 
t Reference  2 for  a qualitative  description  of  the  laminar  flow  recovery), 

the  modeling  employed  in  the  present  report  is  quite  simplified  in 
order  to  accentuate  the  basic  features  associated  with  the  thermal 
recovery  by  a vehicle  after  a loss  of  laminar  flow.  Thus,  instead  of 

* a shaped  body,  we  have  computed  the  recovery  process  on  a flat-plate 
so  as  to  avoid  the  complications  introduced  by  pressure  gradients,  es- 
pecially the  influence  of  pressure  gradients  on  estimating  where  bound- 
ary layer  transition  would  occur  during  the  recovery  stage.  The  flat- 

* plate  configuration  also  allows  fairly  accurate  estimates  of  even  the 
local  transient  heat  transfer  rates  during  the  recovery  process  for 
either  laminar  or  turbulent  flows. 

Another  simplification  which  has  been  employed  herein  is  the  assump- 
tion that  the  vehicle  is  heated  by  an  internal  fluid  flow  system  in 
which  heated  water  (possibly  from  the  vehicle's  propulsion  plant)  is 
pumped  up  to  the  bow  of  the  vehicle  and  then  toward  the  aft  and  in 
contact  with  the  outer  surface.  This  heating  configuration  is  en- 
visioned as  a likely  candidate  from  an  engineering  standpoint  althouqh 
other  heating  schemes  may  be  as  feasible  (viz,  electrical).  Again,  to 
avoid  unnecessary  complications,  the  heating  fluid  passage  is  assumed 


to  be  of  uniform  height  here  despite  the  fact  that  present-day  exper- 
imental vehicles  have  a varying  channel  height  in  order  to  "tailor" 
the  local  heat  flux  along  the  vehicle  surface  to  match  some  other  re- 
aui rements . 

A third  simplification  which  has  been  introduced  into  the  present  an- 
alytical model  is  the  assumption  that  the  vehicle  speed  remains  con- 
stant throughout  the  period  of  laminar  flow  recovery.  Thus,  the  pre- 
sent results  are  applicable  to  the  situation  where  a vehicle  is  towed 
through  the  water  at  a constant  speed.  The  drag  increase  which  accom- 
panies the  loss  of  laminar  flow  and  which  would  cause  a speed  decrease 
in  an  operational  situation  could  be  incorporated  at  some  later  time. 

Such  a consideration  obviously  would  involve  a knowledge  of  the  .ehicle's 
propulsion  system  characteristics  as  well  as  fairly  accurate  estimates 
of  the  parasitic  drag  increments  during  and  subsequent  to  the  loss  of 
laminar  flow.  The  slowing  down  of  the  vehicle  caused  by  a loss  of  lam- 
inar flow  would  have  the  effect  of  shortening  the  laminar  flow  recovery 
time  since  the  natural  boundary-layer  transition  location  and  the  heat 
flux  from  the  vehicle  both  change  with  vehicle  speed.  Thus,  the  present 
estimates  for  laminar  flow  recovery  are  probably  conservative. 

In  order  to  maintain  some  semblance  with  presen.t-day  heated  vehicle 
designs,  a composite  material  outer  wall  for  the  vehicle  has  been  re- 
tained in  the  present  analytical  model.  A uniform  metallic  vehicle  skin 
would  have  been  simpler  to  model  but,  aside  from  the  problems  of  physi- 
cally maintaining  the  integrity  of  such  a skin  (i.e.,  problems  associated 
with  corrosion,  pitting,  wetting,  etc.),  a composite  material  surface 
(such  as  a ceramic  coating  on  a copper  substrate)  has  the  advantage  of 
having  a high  heat  capacity  combined  with  an  insulating  effect  to 
ameliorate  the  chilling  which  accompanies  a sudden  burst  of  turbulent 
flow. 
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Even  with  the  above  simplifications  one  should  appreciate  the  fact 
that  the  transient  heat  transfer  problem  between  the  heating  fluid 
and  the  vehicle  skin  and  between  the  vehicle  skin  and  the  outer  flow 
are  strongly  coupled  to  one  another.  The  derivation  of  the  partial 
differential  equations  which  describe  this  coupling  is  given  in  Section 
2 of  this  report  along  with  the  various  approximations  which  make  the 
solution  of  these  equations  tractable.  In  Section  3,  these  equations 
are  expressed  in  finite  difference  form  so  that  the  wall  and  heating 
fluid  temperatures  can  be  computed  during  the  recovery  period  using 
numerical  techniques. 

Results  from  the  numerical  calculations  are  presented  ir,  Section  4 for 
several  different  scenarios  whereby  laminar  flow  may  be  lost.  These 
scenarios  include:  the  encounter  with  a turbulent  "patch"  embedded  in 
the  ambient  fluid;  the  "sticking"  of  a particle  originally  suspended 
in  the  free  stream  to  the  vehicle  surface  for  a short  period  of  time 
(thus  "tripping"  the  boundary  layer);  and  the  entry  of  a particle  into 
» the  boundary  layer  and  causing  a localized  turbulent  "burst"  which  is 

then  advected  along  the  vehicle  surface.  The  vehicle  parameters  select- 
ed in  order  to  complete  the  numerical  calculations  (i.e.,  free-stream 
speed,  body  length,  skin  thickness  and  composition,  and  heating  water 
flow  rate  and  initial  overheat)  were  based  on  the  characteristics  of 
a class  of  present-aay  tow-tank  models.  Since  no  particular  vehicle 
is  exactly  modeled  and  because  of  the  simplifications  inherent  with  the 
present  modeling,  a one-to-one  comparison  between  the  measured  gross 
operating  conditions  and  those  computed  from  the  present  model  (such 
as  the  total  heat  transferred  or  the  transition  point  location)  would 
probably  not  be  meaningful.  The  laminar  flow  recovery  times  should, 
however,  be  in  the  same  neighborhood  as  might  have  been  observed. 
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A discussion  of  the  numerical  results  (for  relatively  short 
times  into  the  recovery  process)  is  also  contained  in  Section  4.  This 
is  followed  by  an  approximate,  long-time  analytical  solution  for  the 
laminar  flow  recovery.  This  latter  solution,  when  combined  with  the 
numerical  results  allows  an  estimation  of  the  eventual  recovery  time. 

Finally,  conclusions  and  recommendations  for  future  study  are  presented 
in  Section  5 based  on  the  findings  contained  herein. 
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2.  FORMULATION  OF  COUPLED  HEAT  TRANSFER  PROBLEM 

As  indicated  earler,  the  heated  vehicle  is  modeled  as  a flat  plate 
configuration  as  shown  in  the  following  sketch: 


The  vehicle  is  assumed  to  travel  at  speed  U^  in  water  whose  temperature 
is  Heating  fluid  (also  assumed  to  be  water)  is  introduced  into 

the  inner  passage  at  the  leading  edge  at  temperature  T.j  . The  inner 
shell  is  assumed  to  be  an  insulator.  Heat  is  transferred  from  the 
fluid  in  the  inner  passage  to  the  vehicle  skin 'at  a rate  q^(x,t)  and 
is  removed  from  the  vehicle  skin  by  the  external  flow  at  a rate 
q0(x,t).  The  configuration  is,  in  effect,  a co-flow  heat  exchanger. 

The  vehicle  skin  locally  undergoes  transient  heating  or  cooling  ac- 
cording to  whether  q-  is  greater  than  or  less  than  qQ.  The  detailed 
formulation  of  the  analysis  thus  involves  the  coupled  consideration 
of  the  inner  flow  passage,  the  vehicle  skin  and  the  external  flow. 
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2.1  Inner  Passage 

The  inner  passage  is  assumed  to  be  of  constant  cross-sectional  height 
so  that  the  mass  flow  per  unit  area  is  also  constant.  The  flow  in  the 
inner  passage  is  assumed  to  be  a fully-developed  turbulent  flow  of  a 
constant  property  fluid.  All  extensive  quantities  associated  with  the 
heating  water  are  per  unit  span. 

Continuity  considerations  dictate  that  the  mass  flow  per  unit  span, 
m=pua,  is  constant.  However,  sincethe  passage  height,  a,  and  fluid 
density,  p,  are  each  constant,  the  flow  speed,  u,  is  also  constant. 

The  momentum  equation  simply  allows  calculation  of  the  pressure  drop 
in  the  inner  passage  and  is  not  of  direct  concern  to  the  coupled  heat 
transfer  problem. 


The  equation  for  the  fluid  temperature  in  the  inner  passage,  T^(x,t), 
is 

3T.(x,t)  3T.(x,t) 

pac.  + puac.  = - q.(x,t) 

1 3t  1 3x  1 


where  c^  is  the  specific  heat  of  the  heating  water  and  q.  is  the  heat 
flux  from  the  heating  fluid  to  the  vehicle  outer  surface.  Using  the 
traditional  approximation  for  the  heat  flux,  this  equation  is  written 
as : 


9T . 3T . 

— — + u — - 
3 1 8x 


T.j  (x  t.) 


Tw.(x,t) 

l 


(2.1) 


where  Tw.(x,t)  is  the  temperature  of  the  inner  surface  of  the  vehicle 
skin  and  h^  is  the  turbulent  heat  transfer  coefficient  for  the  flow  in 
the  inner  passage.  The  heat  transfer  coefficient  h^  can  be  obtained 
using  Reynolds'  analogy  (see  Reference  3); 


L 


-8- 


hi  = puci  Pr~  /a  , (2.2) 

where  the  friction  factor  f is  a function  of  the  Reynolds  number  based 
on  the  hydraulic  diameter  of  the  inner  passage  and  is  evaluated  from 
information  on  flow  in  circular  pipes  and  Pr  is  the  Prandtl  number  of 
the  heating  fluid. 


» 


* 


l 


Vehicle  Skin 


The  relevant  equation  for  the  vehicle  skin  temperature,  Tw(x,t),  is  the 
transient  energy  equation 


Me 


3Tw(x,t) 

at 


q.(x,t)  - qQ(x,t) 


(2.3) 


In  this  equation,  Me  is  the  heat  capacity  of  the  vehicle  skin  per  unit 
area,  Tw  is  that  skin  temperature  which  best  represents  the  heat  con- 
tent of  the  skin  and  is  described  in  more  detail  below,  and  qQ  is  the 
heat  flux  to  the  outer  flow  from  the  wall.  Because  the  vehicle  skin 
thickness  is  small  compared  to  length  scales  in  the  streamwise  direc- 
tion, we  have  ignored  heat  conduction  along  the  vehicle  in  comparison 
to  conduction  through  the  skin. 


Consider  the  skin  to  be  a composite  of  a metal  with  a thin  outer  ceramic 
cladding,  as  shown  in  the  sketch  below.  The  rationale  for  selecting  a com- 
posite skin  is  described  in  the  Introduction.  The  overall  thickness  is  i. 


The  heat  capacity  per  unit  area  is  then 
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r 


Me  = (pc£)  . , + (pc£)  . (2.4) 

w ' 'metal  v 'ceramic  v ' 

where  here  p is  the  density  of  the  material  under  consideration  and  c is 
its  specific  heat.  For  such  a composite,  the  thermal  resistance,  k,  is 
the  sum  of  the  thermal  resistances  of  each  of  the  two  layers. 


The  effective  thermal  conductivity  of  the  composite  is 


eff 


(t)  ♦(■ 

|)  ( 

"f)  total 

^ 'metal  ' 

' 'ceramic 

(2.5) 


If  the  base  metal  is,  say,  copper  then  the  metal  will  dominate  the  heat 
capacity  of  the  composite,  while  the  ceramic  coating  would  provide  most 
of  the  resistance  to  heat  transfer. 

For  computational  purposes,  it  is  easier  to  work  with  a skin  of  thick- 
ness !i  having  a constant  k ^ as  given  by  equation  (2.5).  In  such  an 
equivalent  skin,  the  effective  metal -ceramic  interface  is  located  ac- 
cording to  their  respective  thermal  resistances  as  shown  in  the  sketch 
below: 


L 


p 


» 


I 
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In  a steady-state  condition  (q^  = q ) the  temperature  decreases  linearly 

in  the  stretched  coordinate  system  as  indicated  in  the  sketch.  For 

heat  capacity  purposes,  the  temperature  T should  be  evaluated  at  or 

w 

near  the  center  of  the  metal  layer,  y = ££.,  where 


£ = 


1 

2 


metal 

total 


(2.6) 


In  a transient  situation  (q^q  ),  the  temperature  gradient  is  not  con- 
stant even  in  a material  of  constant  k As  a first  approximation, 
the  temperature  distribution  through  the  skin  will  be  assumed  to  be 
quadratic  in  the  y-stretched  coordinate  system  as  indicated  in  the 
previous  sketch. 


The  transient  heat  transfer  rates  into  and  out  of  the  vehicle  skin  can 
be  computed  from  Fourier's  law  of  heat  conduction  as: 

V*-1’  * - keff  (E)  <2-7> 

y=0 

qo(x’t)  ‘ - keff  (i)y=t  (2-8) 

where  T(x,y,t)  is  the  local,  instantaneous  temperature  distribution  in 
the  equivalent  skin  with  constant  effective  thermal  conductivity. 


2.2  Outer  Flow 


The  heat  removal  rate 


q0(x,t ) 


can  also  be  written  as 


q0(x,t) 


(2.9) 


where  hQ  is  assumed  to  be  the  local  flat-plate  heat  transfer  coeffi- 
cient. For  a laminar  boundary  layer  (see  Reference  2) 


while  for  a turbulent  boundary  layer  (Reference  3) 


The  boundary  layer  flow  will  be  taken  as  laminar  upstream  of  the  tran- 
sition point  and  turbulent  downstream  of  the  transition  point.  Except 
where  turbulent  flow  is  imposed  (say  by  encountering  a turbulent  patch 
in  the  freestream)  the  transition  Reynolds  number  for  the  flat-plate 
boundary  layer  will  be  assumed  to  be  a function  of  (Tw  - T^)  only  and 

the  calculations  by  Wazzan  et  al.  (Reference  4)  which  are  based  on  the 
g 

e boundary-layer  transition  criterion  will  be  used  to  approximate  that 
functional  dependence. 


3. 


NUMERICAL  METHOD  FOR  SOLVING  THE  COUPLED  EQUATIONS 


As  developed  in  Section  2,  the  transient  coupled  equations  for  the 
heating  fluid  and  vehicle  skin  are: 


3Ti  3T.  ^ qi(x,t) 
3t  3x  pac- 

°TW  _ qi (x,t)  - q0(x,t) 

' t Me 

w 

where  the  heat  fluxes  are 


(3.1) 


(3.2) 


qi 


(3.3) 


q = h 
o o 


Two  - T--’ 


(3.4) 


and  the  transition  location  between  laminar  and  turbulent  flow  states 
is  determined  from  the  approximation  shown  in  Figure  3.1  (page  17). 


The  temperature  distribution  across  the  vehicle  skin  is  assumed  Quadratic 
at  any  instant  and  for  all  axial  locations  x so  that 

T(x»y»t)  = aQ(x,t)  + a j(x,t) y + a2(x,t)y2  (3.5) 


with  boundary  conditions 

T(x,0,t)  = Tw.(x,t) 

r(XjLjt)  ~ tWq ( x , t ) 


(3.6) 


The  coefficients  a0>  a^  and  a^  can  be  obtained  as  follows,  since 
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a = T 
o wi 


a,  = 


a„  = 


£2T  - (£2-l)  T . - T 

wo  - ' wi  w 

; : 

T - ;T  + T Af-  1) 
w wo  wi v ' 

: :2(-'-l) 


(3.7) 


ho  expressions  for  the  neat  f :ux  into  and  out  of  the  skin  are  therefore 


q.  = - k ,,  — , 

eff  svi  n 
- l v=0 


eff  'yi 
^ ( y=' 


eff 


1 eff 


r2T  - ( r 2- 1 )T  . - T 
• wo  - ‘ wi  w 

£ (C-l) 


'(r-2)Two-  (c>1)\i  + T, 

5 (5-1) 


(3.8) 


(3.9) 


•_'Se  e * ' '.ons  can  t • rewritten  as 


■ . = . T + T + • T 
1 ' wo  2 wi  ' 3 ' w 


T -*•  T -*■  T 

1 wo  2 wi  '2  in 


(3.10) 

(3.11) 


where 


eff 


( -1) 


eff  (~-2) 
(C-l) 


o _ ''eff  ( 2+1 ) 


. keff  (c-l) 

’ x7 

i c 


_ eff  1 


(C-l) 


eff  1 


i 5(5-1) 


Note  that 


*1  * 2 + = 0 and  °1  + p2  * s3  = 0- 


-14- 

• 

Let  AT  = T - T , AT.  = T.  - T , AT  = T -T  and  AT  . = 
w w »’  i i oc’  wo  wo  “ wi 

Then,  by  equating  equations  (3.10)  and  (3.11)  with  eouations 

T . - T . 

WI  « 

(3.3)  and 

(3.4),  we  have 

B,  AT  n + (*+h.)  AT  . = h.  AT.  - B,  AT 

1 wo  ' 2 v wi  it  3 w 

(3.12) 

(a.-h  ) AT  + a0  AT  . = - a,  AT 
' 1 o wo  2 wi  3 w 

(3.13) 

The  solution  of  equations  (3.12)  and  (3.13)  for  Tw-j  and  Tw0 

is 

1 

• 

- _ ''  1 " hc  ;hi-'-  ‘ i + /3S 1 " S3^1  " ho\  1Tw 

'wi 

1 - h )l  c + h . ) . a 

1 o " 2 V -2-1 

= g.  ■'  T.  + q„  AT 

31  1 -2  w 

(3.14) 

• 

- _ ”a2hi“ 1 i + _a2£3  ' a3^2  + hi  \ ~Tw 

(al  " ho^~2  + hi ' " a2;"  1 

• 

9o  AT.  + g.  AT 

3 3 i 34  w 

(3.15) 

where 

1 

9j  = (“1  - V VE 

(3.16) 

g2  " [a3£l  ‘ £3^al  " ho^  / E 

(3.17) 

> 

g3  " ‘ a2hi/E 

(3.18) 

d 

I 


I 


9^  - 3 * a3 ( S-2  ^ ^ £ 


(3.19) 
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and 

E = (a1  - hQ)(fc2+hi)  - a261  . (3.20) 

The  coupled  partial  differential  equations  (1)  and  (2)  for  the  depen- 
dent quantities  T_.  and  T can  be  then  rewritten  as 


» 


> 


( 


I 


» 


I 


9AT.  9 (AT.) 

]_  £ 

3t  L 3 (x/L) 


hi 


pac^ 


P-3l)  -T,  - 92  tT„] 


^.(l-sp-V: 


Me 


h .g0  + h g„  AT 
L i 2 o34J  w 


Me 


(3.21) 


(3.22) 


where  L is  the  overall  length  of  the  vehicle. 


The  heat  transfer  coefficients  hQ  and  h^  are  taken  as 


and 


h 


o 


for  laminar  flow 


for  turbulent  flow 


(2.23) 


h. 


rue  . 

l 


fully  turbulent  channel  flow 

(3.24) 


where  the  friction  factor  f is  a function  of  the  Reynolds  number  based 
on  the  hydraulic  diameter  of  the  inner  passage  and  is  evaluated  from  in- 
formation on  flow  in  circular  pipes. 


> 


The  transition  Reynolds  number  versus  the  required  (constant)  wall  tem- 
perature overheat  for  a flat  plate  can  be  plotted  as  shown  in  Figure  3.1. 
For  computational  purposes,  this  curve  is  fit  approximately  by  the  fifth- 
degree  polynomial 


(AT) 


tr 


= c + c, 
o 1 


Re 


tr 


10" 


+ c. 


. .c, 


Re 


tr 


ict 


(3.25) 


where 


c = - 4.31002 
o 

Cj  = 1.29390 

c2  = - 2.78416  x 10' 

c,  = 3.792A9  x 10' 

c4  = - 2.52094  x 10 

c n = 6.91835  x 10 

□ 


By  comparing  the  wall  overheat  temperature  AT  determined  from  equation 
(3.15)  (using  solutions  for  AT  ■ and  AT  ) , and  the  transition  criterion  of 
equation  (3.25),  one  can  determine  the  regimes  of  laminar  or  turbulent  flow. 
Note  that  this  procedure  ignores  the  influence  of  the  upstream  wall  temper- 
ature variations  since  comparisons  are  made  locally  with  a transition  cri- 
terion based  on  an  assumed  uniform  and  steady  wall  temperature.  To  in- 
corporate the  spatial  arid  temporal  changes  of  the  upstream  wall  tempera- 
ture into  a transition  criterion  would  be  so  complicated  and  time  consum- 
ing that  it  would  overwhelm  all  other  considerations.  In  effect,  a 


Wall  lemperaturo  Overheat 


35  - 


: 20- 


. Laminar  Flow 


.Turbulent  Flow 


c = - 4.31002 
0 

c.  = + 1.29390 

c2  = - 2.78416  x 10' 

c0  = + 3.70249  x 10' 


c4  = - 2.52094  x 10'u 
cc  = + 6.91835  x 10'5 


«...  = c. 


_C:  (w) 


Retr\ 


Figu re  3 . 
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complete  boundary-layer  stability  calculation  would  be  required  at 
each  time  step  during  the  recovery  period.  It  was  the  authors'  be- 
lief that  the  above  simple  criterion  would  provid"  reasonably  accur- 
ate estimates  of  the  transition  point  throughout  the  recovery  period 
(at  least  to  within  the  order  of  accuracy  of  other  assumptions  made). 


l 


i 


3. 1 Steady  Flow  Conditions 

When  a stead-state  condition  is  achieved,  equations  (3.21)  and  (3.22) 
may  be  combined  into  a single  ordinary  differential  equation  whose 
closed-form  solution  can  be  derived  once  the  heat  transfer  coefficients 
h^  and  hQ  are  specified.  For  computational  reasons,  however,  it  was 
more  expeditious  to  formulate  the  separate  equations  in  finite-difference 
format  and  solve  the  equations  numerically.  This  allowed  initial  con- 
ditions for  the  transient  solutions  to  be  calculated  directly  (and  with 
similar  algorithms)  rather  than  having  to  be  computed  from  the  (compli- 
cated) closed-form  solution.  It  should  be  added  that  the  numerical  and 
analytical  solutions  agreed  with  one  another  quite  closely.  The  steady- 
state  form  of  equations  3.21)  and  (3.22)  are: 


- ( ' \ 
i J 

‘i  L 

5(x/L) 

u 

pac . 

AT  = 
w 

h • ( 1-g, ) - h g, 
i 1 o33 

hi92+h094 

g,AT 
3 2 w 


(3.26) 


(3.27) 


Equations  (3.26)  and  (3.27),  when  written  in  finite  difference  form, 
become 


h . L 


(AT.)  = (.  T.)  - A(x/L)  — 

j + l J 


ac . u 


(l-g^UV  -g  (AT  ) 

J J-1 


(3.28) 


Mi-aJ  - hn9? 

(AT  ) = — ± (AT.) 

w • . ' i ' . 

J hi92+nog4  j 


(3.29) 


The  initial  condition  for  AT . is  prescribed  and  the  initial  condition 
for  AT  is  competed  from  the  two  companion  relations 


AT  . 
wi  s 


hi  ATi 


x=0  (&,+  h.)  + B,(l-£) 


(3.30) 


at 


wo 


I x=0 


0 


in  conjunction  with  equations  (3.14)  and  (3.15). 


(3.31) 


3.2  Transient  Conditions 


For  transient  conditions,  the  finite  difference  equations  for  the 
heating  fluid  temperature  and  mean  wall  temperature  are 


( ■ t >n+l  _ / T \ ^ 
' “ i j i ; j 


' t 


u 

2A(x/L) 


n 

j + 1 


n 

j-1 


h. 

l 


■ac-j  L 


(l-9l)UTi)3 


(3.32) 


.7  l"*1  = (AT  )n  * At 
W.J  w'j 


hi(1~9l)-ho93 

Me 


- 


h .o0  + h a. 

1-4  q-4  , . .j.  , n 

i. * “ ' w j 


where  the  superscript  n denotes  a value  at  the  nth  time  step  in  the  cal- 
culation and  At  is  the  time  interval. 

The  initial  conditions  for  these  equations  are  identical  to  those  mentioned 
in  Section  3.1  at  x = 0 and  for  x>0  they  are  taken  as  the  steady-state 
solution  described  in  Section  3.1  (usually  the  laminar  flow  solution  ex- 
cept for  the  casp  of  an  encounter  with  an  infinitely  long  turbulent  patch 
in  which  case  the  ful ly-turbulent , steady  state  temperature  distribution 
was  used). 


r 


I 
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3.3  Computer  Program 

A computer  program  based  on  an  explicit  finite-difference  method  was 
written  to  solve  the  coupled  difference  equations.  The  procedure  used 
for  each  calculation  is  described  below: 

1.  Input  the  data  including  initial  flow  state  (laminar  or 
turbulent) . 

2.  Calculate  the  steady  state  temperature  distributions  from 
equations  (3. 28)  and  (3.29). 

3.  Calculate  the  transient  temperatures  by  using  the  finite 
differenced  equations  (3.32)  and  (3.33). 

4.  Use  the  transition  criterion  (equation  (3.25))  to  determine 
the  local  flow  condition  (i.e.,  either  laminar  or  turbulent 
external  flow). 


I 
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4.  EXAMPLES  AND  DISCUSSION 

In  this  study,  four  types  of  transient  response  scenarios  are  treated. 
They  are: 

1.  Recovery  from  a long  turbulent  patch.  In  this  case,  the 
vehicle  is  assumed  to  have  been  exposed  to  fully  turbulent 
flow  for  a long  enough  period  to  reach  turbulent  steady- 
state  conditions. 

§ 

t<_  0 steady  turbulent  external  flow 

t>0  transient  recovery  toward  laminar  external  flow 

2.  Response  to  and  recovery  from  an  encounter  with  a finite 
length  turbulent  patch  embedded  in  the  free-stream. 

t<0  steady  laminar  flow 

-T^ < t < 0 fully  turbulent  external  flow  for  a time  period 
of  Tt  (sec) 

t;  0 transient  recovery  toward ,lami nar  external  flow 


3.  Response  to  and  recovery  from  a temporary  boundary  layer 
trip  at  a fixed  surface  location 

t<^0  steady  laminar  flow 

-Tt<  t<  0\  turbulent  flow  conditions  for  a time  period  of 
an(j  | Tj.  (sec)  downstream  of  (ft) 

| 

t>  0 transient  recovery  toward  laminar  external  flow 
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4.  Recovery  from  a turbulent  spot  generated  at  location 
within  the  laminar  boundary  layer.  This  spot  is 
then  swept  by  the  flow  along  the  flat-plate  surface 
according  to  the  sketch  shown  below. 


A 

( x . + 0 . 5 U t ) 


_ami nat 
FI  ow 


trailing  edge  of  turbulent  patch 

turbulent  /'"'leading  edge  of  turbulent 
flow.  ' patch 

( x t+  U t ) 


/ / 


Laminar  Flow 


These  "turbulent  spots"  are,  of  course,  assumea  to  be 
two-dimensional  in  the  present  context. 


4 . 1 Configuration  Parameters 


The  flat-plate  is  assumed  to  have  a length,  L=  12.45  ft.  The  skin 
composition,  which  is  believed  to  be  representati ve  of  present-day 
experimental  vehicles,  is  a 0.25"  thick  sheet  of  copper  coated  with 
0.007"  of  a ceramic  material.  For  this  skin  composition  and  using 
tabulated  thermal  properties  for  each  material  , such  as  are  given  in 
Reference  5,  it  may  be  shown  that: 


k 


0.495  ■ ■ 

sec-ft2-°F 


Me  = 1.084  -4— 
w ft2  °F 


0.0844 
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The  heating  water  is  assumed  to  have  a flow  speed  of  1.53  ft/sec  in  a 

channel  0.01  ft  thick.  This  corresponds  to  a flow  rate  of  0.96  lb  /sec 

m 

per  foot  of  span  width  which  is  also  believed  representative  of  present 
day  experimental  vehicle  flow  rates. 


i ♦ 

( 

< 

f 
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For  each  choice  of  operating  conditions  (free-stream  speed,  heating 

water  flow  rate  and  inlet  temperature)  the  steady  state  temperature 

distributions  along  the  length  of  the  flat  plate  were  computed.  The 

results  for  a typical  case  (U  = 20  knots,  T_  = 70°F,  (z.T,-)  = 55°F)  are 

0 

shown  in  Figure  4.1.  Mote  that  the  outside  wall  temperatures,  both 
laminar  and  turbulent,  rise  abruptly  from  the  leading  edge  (where  the 
temperature  must  equal  the  free-stream  fluid  temperature)  to  a maxi- 
mum value  at  between  10  and  30.  of  the  flat  plate  length  and  decrease 
slowly  thereafter.  Mote,  however,  the  large  differences  in  level  be- 
tween the  two  cases.  In  order  to  properly  maintain  a steady  laminar 
external  flow,  the  outside  wall  temperature  difference  should  every- 
where exceed  that  prescribed  by  the  e^  transition  criterion. 


In  the  laminar  situation  the  heating  water  loses  44'  of  its  initial 
heat  content  while  in  the  turbulent  case,  because  of  the  higher  heat 
transfer  rate  to  the  external  flow,  the  heating  water  loses  78;  of  its 
initial  heat  content.  At  the  trailing  edge  of  the  plate  with  fully  lami- 
nar flow,  the  outside  wall  temperature  exceeds  that  required  to  just  main- 
tain laminar  flow  by  3.0°F.  If  the  initial  heating  water  temperature  dif- 


at the  trailing  edge  of  the  plate  would  be  reduced  to  zero.  It  will  be 
shown  in  Section  4.4  that  the  time  required  for  laminar  flow  to  be  fully 
recovered  after  some  turbulent  event  is  a sensitive  function  of  this 
temperature  excess.  When  the  wall  temperature  excess  at  the  trailing 
edge  is  zero,  the  recovery  time  is  infinite. 


9 


Temperature  Difference, 


The  results  of  the  steady-state  calculations  for  free-stream  speeds 
of  20,  40  and  50  knots  are  summarized  in  Table  4.1.  The  changes  of 
the  laminar  heat  loss  and  the  trailing  edge  wall  temperature  excess 
to  changes  in  the  initial  heating  water  temperature  (above  the  mini- 
mum required  value)  are  shown  in  Table  4.2. 


Tabl  e 


Table  4.2 


SENSITIVITY  TO  HEATING  WATER  INLET  TEMPERATURE  DIFFERENCE 

(L  = 12.45  ft.  m = 0.96  lb  /sec  - ft  of  width) 

in 


Free-Stream  Speed 

( knots ) 

' 

Laminar  Heat  Flux  Increment 
per  unit  Heating  Water 
Temperature  Increase 

Trailing  Edge  Wall  Temp. 
Excess  per  unit  Heating 
Water  Temp.  Increase 

kw/ft  Of  width  per  °F 

°F/°F 

1 

20 

1 

0.47 

I 

0.43 

40 

0.53 

0.34 

50 

0.53 

0.32 

4 . 3 Results  o_f_  Numerical  Calculations 

The  results  of  the  calculations  for  the  various  recovery  scenarios  will 
now  be  presented.  The  results  will  be  shown  only  for  about  the  first  four 
seconds  into  the  recovery  history.  The  principa-1  reason  for  this  is  that, 
for  later  times,  the  explicit  finite-difference  technique  used  to  solve 
the  coupled  partial  differential  equations  becomes  unstable.  This  numer- 
ical instability  manifests  itself  initially  as  small  "wiggles"  about  the 
mean  axial  distribution  after  about  8 seconds  into  the  recovery  period 
and  these  wiggles  become  larger  and  larger  at  later  times.  This 
instability  is  a result  of  the  numerical  algorithm  employed  rather  than 
being  a result  of  the  equations.  Rather  than  attempting  to  incorporate 
a more  sophisticated  technique  in  order  to  extend  the  computations,  we 
have  instead  found  an  approximate  analytical  solution  which  is  valid  for 
the  latter  stages  of  the  recovery  period.  This  analytical  solution  is 
described  in  Section  4.4  of  the  present  report.  The  key  to  combining 


the  early-time  numerical  solution  with  the  late-time  analytical  solu- 
tion is  that  the  latter  requires  an  initial  condition  which  is  provided 
by  the  numerics  at  four  seconds.  Thus,  the  entire  recovery  process  is 
describable  and  the  matching  between  the  early  and  late  time  solutions 
is  not  critically  dependent  on  the  precise  matching  time.  The  details 
of  the  matching  procedure  and  estimates  for  the  recovery  time  are  also 
described  in  Section  4.4. 


4.3.1.  Recovery  From  "Long"  Turbulent  Patch 

The  initial  wall  and  fluid  temperatures  used  for  the  recovery  calcula- 
tions described  in  this  section  are  taken  as  those  corresponding  to  a 
ful ly-turbul ent  external  boundary  layer.  This  is  the  most  "severe"  de- 
viation imaginable  from  a fully  laminar  external  boundary  layer  and 
corresponds  to  the  flat  plate  encountering  a very  long  turbulent  patch 
embedded  in  the  ambient  fluid.  Three  representative  calculations  are 
shown  in  Figures  4.2,  4.3  and  4.4  where  the  outside  wall  temperature 
histories  are  shown  for  several  seconds  after  the  flat-plate  exits  from 
the  turbulent  patch. 

In  Figure  4.2,  for  example,  the  steady-state  laminar  and  turbulent  flow 
wall  temperature  distributions  are  shown  as  well  as  the  wall  temperature 
distribution  required  to  prevent  natural  boundary  layer  transition  accord- 

9 

ing  to  the  e'  criterion.  The  outside  wall  temperature  distributions, 

' Tw  , for  several  times  after  exit  from  a turbulent  patch  are  shown.* 

• * r\ 


* The  present  calculations  assume  that  at  t=0  the  imposed  external  tur- 
bulent flow  completely  vanishes  leaving  only  the  natural  transition  to 
turbulent  flow  within  the  boundary  layer  rather  than  accounting  for  the 
fact  that  the  turbulent  patch's  edge  would  move  rearward  along  the  plate 
at  33.8  ft/sec.  The  difference  between  these  two  is  slight  since  the 
natural  transition  point  moves  so  much  slower  than  the  free-stream 
speed  (this  was  the  case  in  every  calculation). 
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Note  that  after  removal  of  the  imposed  external  turbulent  flow  in  Figure 

4.2,  the  transition  location  virtually  "jumps"  to  a location  30%  aft  of 

► the  leading  edge  and  then  progresses  at  a much  slower  rate  toward  the 

r 

trailing  edge  (i.e.,  in  the  first  0.05  seconds,  the  transition  location 
moves  from  x/L=0.2  to  x/L~0.3  and  nearly  four  seconds  later  the  tran- 
sition point  has  only  progressed  to  x/L~0.4).  This  characteristic  “jump 
and  creep"  motion  of  the  transition  point  has  been  observed  in  nearly 
every  calculation.  The  heating  water  temperatures  are  also  shown  in 
Figure  4.2  (only)  for  the  fully  laminar,  fully  turbulent  and  for  an  inter- 
mediate time  during  recovery. 

i 4.3.2  Recovery  From  Finite  Length  Turbulent  Patch  Encounter 

If  the  flat-plate  were  to  encounter  a finite  length  turbulent  patch  em- 
bedded in  the  ambient  water, the  wall  temperature,  upon  exiting  the  patch, 

, would  not  be  as  depressed  as  in  the  previous  case.  Because  of  this,  one 

might  expect  that  the  recovery  time  would  be  shorter  than  the  case  con- 
sidered in  Section  4.3.1.'  To  demonstrate  that  this  is  indeed  the  situ- 
ation, calculations  for  operating  conditions  similar  to  those  considered 

• previously  were  performed  but  for  which  the  turbulent  "soak  time"  was 
several  seconds.  The  calculated  outer  wall  temperature  histories  for 
two  of  these  cases  are  shown  in  Figures  4.5  and  4.6.  For  both  of  these 
the  turbulent  patch  length  was  taken  as  67.6  feet  (i.e.,  a duration  of 

• 2 seconds  and  1 second  for  speeds  of  20  knots  and  40  knots,  respectively). 
After  exiting  the  turbulent  patch,  the  natural  boundary  layer  transition 
point  is  again  seen  to  "jump"  to  an  initial  location  and  then  to  "creep" 
toward  the  flat-plate's  trailing  edge.  Notice  that  downstream  of  the 
transition  location,  the  wall  temperature  continues  to  decrease  even 
after  the  turbulent  patch  is  exited  because  the  boundary  layer  is  still 
turbi  ent  there,  and  as  long  as  it  remains  turbulent,  it  tends  toward 
the  level  for  steady  turbulent  flow.  The  eventual  laminar  flow  recov- 
ery time  for  these  (and  other)  cases  will  be  discussed  in  Section  4.4. 


Figure  4.5.  Recovery  From  Encounter  with  Turbulent  Patch 
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Figure  4.6.  Recovery  From  Encounter  with  Turbulent  Patch 
U = 40  Knots  (67.6  ft/sec) 
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4.3.3  Recovery  From  a Particle  Temporarily  "Stuck"  to  the 
Vehicle  Surface 


One  other  situation  whereby  laminar  flow  could  be  lost  is  by  the  ad- 
herence of  a solid  particle  (say,  originally  suspended  in  the  ambient 
fluid)  to  the  surface  of  the  heated  body  which  then  "trips"  the  laminar 
boundary  layer.  This  situation  was  also  numerically  modeled  and  some 
representative  results  are  shown  in  Figures  4.7  and  4.8  for  the  same 
two  operating  conditions  considered  in  the  previous  section.  For  these 
figures,  it  has  been  assumed  that  the  particle  adheres  to  the  flat-plate 
surface  two  feet  aft  of  the  leading  edge  and  that  it  remains  stuck  for 
a vehicle  travel  distance  of  20.3  feet  (0.6  sec.  duration  at  20  knots, 

0.3  sec.  at  40  knots).  Because  this  particular  case  is  less  severe 
than  the  previous  two  cases,  the  natural  transition  point  almost  "jumps" 
past  the  plate's  trailing  edge  once  the  particle  is  dislodged  as  may  be 
seen  in  both  Figures  4.7  and  4.8.  After  this  initial  jump,  the  previously 
described  "creep"  toward  the  trailing  edge  is  also  noted,  however. 


4.3.4  Recovery  From  a Turbulent  Spot  Created  Within  the  Original 
L aminar  Boundary  Layer 

The  last  laminar  flow  loss/recovery  history  computation  is  representa- 
tive of  the  situation  where  a particle  (originally  suspended  in  the  am- 
bient fluid)  enters  the  flat  plate's  laminar  boundary  layer  and,  in  doing 
so,  locally  disrupts  the  flow  field.  In  effect,  a turbulent  spot  is 
created  which  is  then  swept  by  the  flow  along  the  length  of  the  plate. 

The  leading  edge  of  the  turbulent  spot  is  assumed  to  be  advected  at  the 
free-stream  speed  whereas  its  trailing  edge  is  assumed  to  move  down- 
stream at  one-half  the  freestream  speed.  Thus,  the  spot's  extent  grows 
as  it  is  swept  along  the  plate. 

Calculations  of  such  an  occurrence  are  shown  in  Figures  4.9  and  4.10 
again  for  the  two  nominal  operating  conditions  considered.  In  either 
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of  these  figures,  the  motion  and  streamwise  growth  of  the  turbulent 
spot  can  be  easily  visualized  as  can  the  effect  on  the  wall  temperature 
(the  step-like  appearance  of  the  computed  wall  temperature  distribution 
is  not  a manifestation  of  the  numerical  instability  mentioned  previously 
but,  instead,  is  an  artifact  of  the  step-like  manner  - both  in  time  and 
position -by  which  the  finite  difference  equations  are  solved).  Note 
that  in  both  of  the  cases  shown  in  Figures  4.9  and  4.10,  the  laminar 
boundary  layer  flow  is  fully  re-established  as  soon  as  the  turbulent 
spot  is  swept  past  the  trailing  edge. 

If  the  particle  entry/turbul ent  spot  generation  mechanism  is,  in  fact, 
one  whereby  laminar  flow  could  be  temporarily  lost  (and  there  is  some 
experimental  evidence  which  substantiates  this  surmis  ) then  it  is  very 
unlikely  that  only  one  such  occurrence  would  be  encountered.  Instead, 
particles  would  be  encountered  at  random  time  intervals  distributed 
about  a mean  which  is  characteristic  of  their  undisturbed  spacing  arid 
the  vehicle  speed.  Thus,  a more  realistic  simulation  of  this  mechanism 
would  be  to  compute  the  effect  of  repetitive  turbulence  spots  being 
swept  over  the  flat  plate  surface.  This  was  not  done  in  the  present 
study  since  a mean  arrival  rate  could  not  be  confidently  estimated  and 
because  the  two  limiting  cases  of  either  a single  particle  or  a contin- 
uous "swarm"  of  particle  arrivals  are  contained  in  the  calculations  al- 
ready described. 


4 . 4 Time  for  Recovery  o f Laminar  Flow 

The  process  of  recovering  a fully  laminar  flow  after  experiencing  tur- 
bulent events  such  as  those  described  in  Section  4.3  can  be  viewed  in 

terms  of  the  approach  of  (LTW)  to  its  steady  laminar  value.  Consider 

0 

the  temperature  distributions  in  the  following  sketch: 
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Temperature  distribution  1 corresponds,  say,  to  that  at  one  of  the  cal- 
culated times.  Temperature  distribution  2 can  be  considered  to  be  that 
at  the  last  time  where  a 'reliable  calculation  is  obtainable.  Tempera- 
ture distribution  4 is  that  for  steady  external  laminar  flow,  and  tem- 
perature distribution  3 is  the  earliest  one  for, which  the  transition 
point  reaches  the  end  of  the  plate  (x/L  = 1.0).  The  recovery  time  can 
be  estimated  from  the  rate  at  which  temperature  profiles  1,  2 and  3 ap- 
proach prof i 1 e 4 . 


At  any  station  x/L  where  the  external  flow  is  laminar  (indicated  by  the 

vertical  dotted  line  in  the  sketch),  equations  (3.21)  and  (3.22)  for 

AT . and  AT  can  be  written 
i w 


9(ATi) 

3t 


9 (AT.) 


3 (AT  ) 
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H1  (AV  + H2  (AV 


(4.1) 
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where 
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H,  = 
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Near  equilibrium  (here  represented  by  profile  4) 
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ATi  = (ATi)  + 5 (t) 
eq 


AT  = (AT  ) + 6 (t) 

w w'eq  wv  ' 


(4.4) 


Substituting  relations  (4.4)  into  equations  (4.1)  and  (4.2)  yields  the 
following  two  sets  of  equations: 

For  the  equilibrium  distribution: 


3(AT.) 

T ^ = - H.  (AT.)  + H (AT  ) 

/-v  1 1 eq  2 w eq 


!x 
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H (AT  ) = H (AT  ) 

j i eq  ^ w eq 


V 
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(4.5) 


For  the  approach  to  equilibrium: 


35. 

T_ 

3t 


H,  6.  + H0  5 
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= h3  fii 


H,  6 
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Equations  (4.5)  are,  in  fact,  the  equations  for  steady  laminar  flow  whose 
solution  would  yield  the  temperature  profile  denoted  as  4 in  the  previous 
sketch.  Since  the  coefficients  H^,  , H.  and  are  all  constant,  we 

seek  solutions  to  equations  (4.6)  that  are  exponential  in  time. 


Let 


vv 


-At 


(4.7) 


Equations  (4.6)  then  become 


(Hf  A)6i  - H26w  - 0 | 

/ (4.8) 

h3  4i  - <H4-X>5w  ' 0 ) • 

For  a non-trivial  solution  to  these  equations  to  exist,  the  quantity  A must 
•satisfy  the  following  equation 

A2  - (H1  + H4)A  + (H1H4-  H2H3)  = 0 • (4.9) 


The  calculated  values  of  A for  the  three  vehicle  speeds  considered  are: 
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Table  4.3.  Characteristic  Time  Constants 


Freestream  Speed 
( knots) 

X1  i 

(sec  *) 

T1 

(sec) 

X2_i 

(sec  *) 

t2 

(sec) 

20 

0.020 

50 

0.547 

1.83 

40 

0.027 

37 

0.555 

1.80 

50 

0.030 

33 

0.558 

1.79 
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Only  and  its  reciprocal,  the  time  constant  Tj , are  relevant  to  the 
present  problem  since  for  • = 0 , the  temperature  deviations  6-  and  6 

are  opposite  in  sign  which  is  not  appropriate  for  the  recovery  processes 
considered  in  this  paper.  For  a = , 6^  and  6 have  the  same  sign  and 

in  fact,  for  all  the  cases  considered  herein,  the  equilibrium  laminar 
wall  and  heating  fluid  temperatures  are  both  approached  from  below. 

A.n  analysis  of  the  root  shows  that  to  a very  close  approximation 

(Me  + oac • ) 

-i  i.i2 S5! — (4.10) 

h 

o 

indicating  that  the  time  constant  is  directly  proDortional  to  the  heat 
capacity  per  unit  area  of  the  skin  and  heating  fluid  together  and  in- 
versely proportional  to  the  heat  loss  rate  to  the  external  flow. 


he  procedure  for  estimating  the  recovery  time  follows  directly  from 
Equation  14.7)  and  the  sketch  at  the  beginning  of  this  section.  One 
may  write 


w 


- (AT  ).  . 

wQ  laminar 


(4.11) 
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More  specifically,  for  the  final  calculated  profile  (profile  2)  and 
the  recovery  profile  (profile  3)  we  may  write  respectively 
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(AT)?  = A e ll 

t 

recovery 

(-T)t>e_  = (£T)3  = 4 e T1 


or  taking  the  ratio  of  the  two  temperature  differences 
(AT). 


'T.E.  _ / 
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^recovery  ~ "2^ 
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(4.12) 


where  (AT)^  = (AT)t  ^ is  the  trailing  edge  wall  temperature  excess. 
Solving  Equation  (4.12)  for  the  recovery  time  yields 


t 


recovery 


= t , In 


'(at)2  ‘ 


+ lr 


(4.13) 


Note 


that  when  ( .IT ) -j- - goes  to  zero,  the  recovery  time  becomes  infinite. 


Recovery  times  have  been  estimated  using  Equation  (4.13)  for  all  the 
cases  computed.  Selected  results  related  to  the  cases  presented  in 
Section  4.3  are  riven  in  the  table  below: 


Table  4.4 


Recovery  from  a Finite  Length  ~urbulent  Patch  Encounter 


Figure  4.2 
Figure  4. 


Figure  4.4 


Figure  4.5 

Figure  4.6 

Figure  4.7 
Figure  4.8 
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Recovery  times  are  not  given  for  the  cases  of  the  passage  of  a turbu- 
lent spot  (Figures  4.9  and  4.10)  since  as  has  been  noted  earlier,  lami- 
nar flow  is  recovered  as  soon  as  the  turbulent  spot  is  swept  past  the 
trail ing  edge. 

In  all  cases,  laminar  flow  is  eventually  recovered  over  the  entire  plate. 
As  expected,  the  recovery  time  is  a sensitive  function  of  the  trailing 
edge  wall  temperature  excess  and,  therefore,  also  of  the  corresponding 
heating  water  inlet  temperature  excess.  Recovery  times  can  be  reduced 
by  providing  a greater  temperature  excess. 

Recovery  times  from  exposure  to  finite  length  patches  of  turbulence 
are,  of  course,  less  than  those  for  infinite  patch  lengths,  but  not 
a great  deal  less.  This  indicates  that  if  a patch  of  turbulence  ex- 
tends for  the  order  of  three  or  more  vehicle  lengths,  significant  surface 
cooling  occurs  because  of  the  turbulent,  rather  than  laminar,  external 
heat  loss  rate. 

From  momentary  trips  or  turbulent  spots,  recovery  is  either  instanta- 
neous, or  else  the  recovery  times  are  fairly  short,  even  shorter  than 
indicated  by  the  present  calculations.  This  is  because  real  trips  and 
spots  are  three-dimensional  (of  finite  spanwise  extent)  rather  than  two- 
dimensional  as  assumed  in  the  calculations.  The  drag  and  heat  transfer 
consequences  of  such  disturbances  are  dependent  more  on  their  residence 
time  (or  flow  time)  and  are  only  weakly  dependent  on  the  characteristic 
recovery  times,  t^,  of  Table  4.3. 
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5 . CONCLUSIONS  AND  RECOMMENDATIONS 

The  maintenance  of  steady  laminar  flow  on  a heated  submerged  vehicle 
requires  that  the  external  skin  temperature  at  every  station  on  the 
vehicle  be  equal  to  or  greater  than  that  required  to  delay  transition 
to  that  station.  The  ability  of  the  vehicle  to  recover  laminar  flow 
after  a turbulent  event  is  sensitive  to  the  amount  by  which  the  wall 
temperature  exceeds  the  minimum  required  for  laminar  flow.  If  there 
is  no  temperature  excess,  the  recovery  time  will  be  infinite.  When 
there  is  excess  temperature,  the  recovery  time  will  depend  on  the  ex- 
tent to  which  the  turbulent  event  has  lowered  skin  temperatures,  and  is 
characterized  by  a time  constant  that  is  directly  proportional  to  the 
heat  capacity  per  unit  area  of  the  skin  and  heating  fluid  together  and 
inversely  proportional  to  the  laminar  heat  loss  rate  to  the  external  flow. 

For  the  calculated  cases,  if  a patch  of  turbulence  extends  for  the  order  of 
three  or  more  vehicle  lengths,  the  recovery  time  for  reasonable  temperature 
excess  is  of  the  order  of  one  to  three  characteristic  times  (1/2  to  2 min- 
utes for  the  cases  calculated).  For  momentary  trips  or  turbulent  spots,  re- 
covery occurs  very  soon  after  the  spot  or  trip  is  swept  past  the  trailing  edge. 

While,  for  simplicity,  the  vehicle  of  the  present  study  has  been  modeled 
as  a two-dimensional  flat  plate,  the  formulation  is  readily  extendable  to 
the  geometries  of  actual  test  vehicles.  It  is  strongly  recommended  that 
such  analyses  for  actual  geometries  be  carried  out  as  part  of  the  thermal 
design  process  of  prospective  vehicles  and  as  an  aid  in  the  interpreta- 
tion of  test  data.  Attempts  should  be  made  to  accommodate  properly  the 
three-dimensional  character  and  azimuthal  distribution  of  turbulent  spots 
and  momentary  boundary  layer  trips. 
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